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AN  APPROXIMATE  THEORY  OF  INELASTIC  COLLISIONS** 

by 

Jn  mt'S  I-'.  Hor  nig  and  J.  O Hirschfelder 

ABSTRACT 

Approximate  trajectories  fot  describing  the  inelastic  encounte r s between 
rigid  diatomic  Inolecules  are  described  These  trajectories,  together  with 
quantum  mecharvcal  perturbation  theory  are  used  to  derive  an  approximate 
formulation  of  the  inelastic,  differential  scattering  cross  sections  for  such 
collisions.  The  quantum  mechanical  correspondence  principle  is  used  when 
necessary  to  bridge  the  gap  between  classical  and  quantum  descriptions  of  the 
problem. 


**This  work  was  carried  out  at  the*  University  of  Wisconsin  Naval  Research 
Laboratory  under  Contract  N7onr- 28511  with  the  Office  of  Naval  Research. 
One  of  us  (JFH)  wishes  to  acknowledge  the  financial  support  which  made  this 
work  possible.  An  Atomic  Energy  Commission  Predoctoral  Fellowship 
administered  through  the  Oak  Ridge  Institute  of  Nuclear  Studies  provided 
support  while  most  of  this  work  and  all  of  the  preliminary  work  was  being 
done.  A General  Electric  Company  Predoctoral  Fellowship  provided 
8’ipaort  while  this  work  was  be-ng  completed. 
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A theoretical  predit  tion  of  any  of  the  non  equilibrium  properties  of 
dilute  gases  depends  on  detailed  information  about  the  dynamics  of  a binary 
collision.  One  can,  for  example,  calculate  the  transport  coefficients  of  a 
rare  gas  if  one  knows  the  angle  of  deflection  in  an  encounter  as  a function  of 
the  initial  energy  and  angula  r momentum  of  the  two  body  system  For  more 
complex  molecules,  the  probability  of  internal  transitions  accompanying  the 
collision  must  also  be  known.  Wang  Chang  and  Uhlenbeck^  have  derived  a 
formal  treatment  of  the  kinetic  theory  of  polyatomic  molecules  which  is 
expressed  in  terms  of  a set  of  such  differential  scattering  cross  sections,  one 
for  every  possible  internal  transition  Similarly,  an  exact  theory  of  the 
pressure  broadening  of  microwave  spectral  lines  will  require  a knowledge  of 
the  probability  of  internal  transitions  accompanying  collisions 

When  molecules  interact  according  to  spherically  symmetric  laws  of 
force  and  undergo  only  elastic  collisions,  the  complete  dynamics  of  a collision 
can  readily  by  obtained.  If  one  or  both  of  these  conditions  are  not  satisfied, 
the  problem  becomes  very  much  more  difficult.  Several  investigations  of 
transport  phenomena  of  special  molecular  models  with  internal  degrees  of 
freedom  have  been  carried  out,  but  the  two -fold  difficulty  of  solving  for  the 
classical  dynamics  of  an  encounter,  then  of  modifying  the  classical  Boltzmann 
equation  has  usually  led  to  inconclusive  results  A brief  summary  of  such  work 
may  be  found  in  chapter  seven  of  The  Molec  ular  Theory  of  Gases  and  Liquids 


1 


Wang-Chang-Chang,  C.  S.  , and  Uhlenbeck  , G E.  Univ.  of  Michigan 
publication  C M - 68  1 (1951) 
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bv  H’  - sr  hfe  ide  r Curtiss,  and  Bird. 

Rigorous  quantum  mechanical  treatments  of  the  problem  are  in  progress 
and  promise  to  \ield  the  most  fruitful  results  The  scattering  of  spherical 
molecules  from  slightly  elongated  ones  has  been  worked  out  by  Gioumousis  and 
Curtiss2  and  is  being  used  as  the  basis  for  a theory  of  the  interaction  of  diatomic 
molecules.  Buckingham  and  Delgarno^  recently  gave  a remarkably  successful 
solution  for  the  diffusion  and  excitation  transfer  of  metastable  helium  in  normal 
gaseous  helium.  Although  m-lastic  this  helium  - metastable  helium  encounter  is 
a case  of  exact  resonance  with  a spherically  symmetric  potential,  so  that  the 
theory  is  considerably  simplified. 

Unfortunately,  the  exact  treatment  of  more  complicated  encounters  may 
not  be  available  for  some  time,  so  it  would  appear  very  desirable  to  develop  an 
approximate  theory  of  such  inelastic  encounters.  We  might  expect  that  such  an 
approximate  theory,  containing  the  essential  features  of  inelastic  collisions, 
would  at  least:  1)  provide  a basis  for  sc  mi -empirical  prediction  of  transport 

coefficients,  and  2)  indicate  which  features  of  the  inelastic  encounter  cause  important 
modifications  in  macroscopic  phe  nomena. 

In  this  report  we  describe  such  an  approximate  theory  and  apply  it  to 
the  theory  of  transport  properties.  For  convenience  we  will  usually  talk  about 
the  interaction  of  diatomic  molecules,  but  wherever  possible,  we  will  proceed 
from  the  most  general  conmde i ations  so  that  the  final  result  will  at  least  suggest 


2 W is.  ORD,  OOR  8. 

Proc.  Roy.  Soc . A,  -1  3.  506  (1952). 
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a treatment  of  more  c omplif  ated  situations  The  purpose  of  the  arguments  will 
be  two-fold:  first  to  establish  the  validity  of  the  overall  procedure  to  justify 

the  use  of  the  correspondence  principle  which  makes  the  procedure  possible, 
and  secondly  to  indicate  approximations  which  will  make  the  method  feasible 
for  practical  calculations.  We  will  attempt  to  identify  clearly  the  nature  of  all 
approximations  as  they  are  made,  and  to  indicate  their  range  of  validity. 

Beginning  with  the  general  quantum  mechanical  formulation  of  the  problem, 
we  will  set  up  a solution  in  terms  of  a perturbation  picture.  It  will  be  argued  that 
the  translational  motion  of  the  molecules  can  be  treated  classically,  so  that  the 
zero  order  solution  of  this  part  of  the  problem  will  simply  be  classical  elastic 
trajectories.  For  convenience,  various  approximations  to  the  classical 
trajectories  will  be  considered  and  compared.  The  nature  of  a trajectory  in  the 
vicinity  of  a sudden  inelastic  transfer  will  be  discussed  in  general  terms,  giving 
certain  restrictions  or  "selection  rules"  on  the  possible  nature  o:  •'  e tr  • 'sf  r 
at  a point  in  space  . Consistent  with  these  selection  rules,  1 mi  t u.  «- 1 1 ; >;• 
developed  for  treating  an  inelastic  trajectory  as  a perturbation- induced  transition 
between  elastic  trajectories.  Finally,  the  formalism  of  a tentative  procedure  for 
determining  the  inelastic  differential  scattering  cross  sections  will  be  developed. 


M AM-  %/*'*'" 
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A.  QUANTUM  MECHANICAL  FORMULATION  AND  SEPARATION 

OF  COORDINATES 

Consider  a system  of  two  diatomic  molecules  "a"  and  "b",  made  up 
o ( atoms  1 and  2 in  a , and  3 and  4 in  b . For  an  arbitrary  fixed  set 
of  coordinates,  the  Schroedinger  equation  of  the  system  is 

-%r(v-E)llf  ~0  (l, 

t -/ 

Here  ^7;V  is  the  Laplacian  in  the  coordinates  of  the  i'th  atom.  If  we 
introduce  relative  coordinates: 


(/*T\ , -t  rrr\  v)  ( /YY\  3 t /T>1  vj 

equation  (1)  becomes; 

T^f  Vr  Iff  + 77  Ijf  + J/Z  'W 

1-  7k  Vi  V - ^ ( v - £)  4T  = o 
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Here  we  have  introduced  reduced  mass  and  total  mass  symbols 


H*- 


/?a«  rn\x 

'On , f-  /m  ^ 


'ryi  -j  /m  y 
^b  t'Oyiy 


ft  _ ( OT},  ( /?7)3  y ^Y)y) 

''yft  , i-  /~rn  ^ *-  sm  , -r  /yt,  y 

— y?n , -h  -t  nrt-i  -t  /Wy 


This  non  singular,  linear  t r anstor  mat  ion  amounts  simply  to  using  the 
a-  1 1 k ruv*  ti  paration  of  variables  for  two  particle  problem  three  times,  once 
on  each  diatomic  molecule,  and  om  <•  on  the  system  of  two  molecules. 

With  the  Schroedinger  equation  in  this  form  it  is  quite  easy  to  sec 
the  effects  of  various  simple  forms  of  the  potential.  Suppose  the  potential 
is  the  sum  of  separate  terms; 


v = V,rs)  . V,rbl  . V,R,„  (4) 

Then  there  is  a solution 

r.r.)  - Xi:  *?(&.)  for-!  &(**-■:  151 

whe  re 

tV  VI  X + X ~ 0 
7k  & ~ & ( V~‘  ' ££}  & 

-fa  vL  #b  - V(-;  " 

9 - tX  - E«»)  ® =0 

With  this  approximation.  then  we  have  separated  the  problem  into  four 
separate  physically  resonable  pa  ts  Fquation  f> . 1 describes  the  motion  of 


the  center  of  mass  of  the  whole  syste  m moving  like  a free  particle  of  mass  M 
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The  second  and  third  equations  describe  internal  motions  oi  the  molecules 
relative  to  their  mass  tenters,  and  the  last  < qa-ition  ciesc  nbes  the  relative 
motion  of  th”  two  molecules  relative  to  their  joint  center  of  mass.  This 
motion  can  be  described  in  terms  of  the  equivalent  one  particle  problem  in 

„ A la  Me 

which  a particle  oi  mass  /y  - ~pj — moves  under  the  influence  ot  a 
potential  V(R  j . We  wil1  use  this  picture  of  the  equivalent  one  particle 
problem  in  the  rest  o'  this  report. 

The  influence  of  asymmetry  in  the  potential  field  as  far  as  modifying 
the  trajectory  and  inducing  internal  transitions  may  now  be  treated  as  a 
perturbation  problem.  If  the  molecules  had  electric  dipole  moments,  we  would 
have  such  a potential  of  considerable  physical  interest  Before  co  sidering  the 
effect  of  a perturbation,  however  let  us  consider  equations  6 in  i ore  detail. 

The  motion  of  the  center  of  mass  of  the  entire  system  will  in  general 
not  enter  into  kinetic  theory  problems  so  only  the  last  three  equations  will 
b«  of  interest.  Equation  o.4  the  relative  motion  of  the  molecules  can 
frequently  be  replaced  ay  its  classic  al  limit.  The  DeBroglie  wave  length  of 
a molecule  with  energy  k F is  given,  at  room  te mpe r atur e , by  — 

VTT  ^ 

where  M is  the  molecular  weight.  Or,  from  the  point  of  view  of  the 
uncertainty  principle  an  uncertainty  of  1 X in  the  position  leads  to  an 


t ^ \r  • 0^3 

uncertainty  in  the  velocity  such  that  — rzr  — . j again  at  room  temperature. 
Thus,  to  a reasonable  approximation,  a classical  trajectory  may  be  used  to 
describe  the  relative  motion  of  the  molecules  This  approximation,  of  course, 
is  the  basis  of  the  entire  classical  theory  of  transport  properties;  its  validity 


over  a rather  wide  range  of  conditions  fs  shown  by  the  success  ol  that  theory. 
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Clearly  no  such  classical  limit  may  be  taken  for  equations  6. 2 and 
6.3.  Solved  with  the  usual  approximations,  these  equations  give  the  harmonic 
oscillator  and  rigid  rotor  wave  functions.  At  ordinary  temperatures  only  the 
lower  vibrational  and  rotational  levels  will  be  populated,  so  that  the  quantum 
mechanical  description  of  these  coordinates  of  the  system  is  necessary.  In  a 
great  many  special  cases  it  will  turn  out  that  excitation  of  vibrational  levels 
is  very  improbable  in  thermal  collisions,  so  that  only  the  rotational  motion 
need  be  considered.  For  this  reason  we  will  concentrate  on  a treatment  of 
the  rotational  motion,  though  the  extension  to  include  vib-ational  transitions 
would  not  be  difficult. 


B.  APPROXIMATE  TRAJECTORIES  FOR  ELASTIC  COLLISIONS 

Since  the  effect  of  a perturbation  will  in  general  be  viewed  as  causing 
transitions  among  various  unperturbed  states  of  the  system,  it  will  first  be 
necessary  to  investigate  these  unperturbed  states.  For  the  internal 
coordinates  of  a diatomic  molecule,  these  are  obviously  the  quantum  mechanical 
wave  functions  for  the  harmonic  oscillator  and  the  rigid  rotor.  Since  in  the 
unperturbed  state  these  do  not  change  during  a collision,  the  remainder  of  the 
zevoth  order  approximation  is  just  the  classical  trajectory  for  the  elastic 
collision  of  spherically  symmetric  molecules. 
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For  the  elastic  classical  trajectory  we  must  first  decide  on  some  sort  of 
effective  spherical  potential  which  will,  in  general,  depend  on  the  internal  states 
of  the  two  molecules.  In  the  case  of  eiectric  dipoles,  this  could  be  done  by 
calculating  the  interaction  energy  of  a pair  of  "clamped"  rotating  dipoles.  This 
calculation  has  been  done  by  London  with  his  "Unscharfe  Resonance",  a modified 
second  order  perturbation  technique.  In  any  case,  we  will  assume  for  the 
mement  that  we  have  such  a potential  V (Ra^)  * ^^ere  1 and  j denote  the 

total  quantum  states  of  molecules  a and  b , and  is  the  inter  molecular 

distance.  We  also  assume  that  this  potential  contains  the  necessary  hard  core 
or  repulsive  properties 

Elastic  trajectories  might  now  be  calculated  in  the  usual  manner. 

Spec ification  of  g , the  initial  relative  velocity,  and  b . the  collision  parameter 

or  miss  distance  gives  the  energy  and  the  absolute  value  ot  the  angular  mementum 

of  the  system.  The  distance  of  closest  approach,  r can  now  be  calculated 

m 

as  the  greatest  positive  root  of  the  equation 


| _ 


r v 


0 


(?) 


Then  the  angle  of  deflection  % is  given  by 

X * 7r  - 2 - rv  - H&to 


. !_ 

X 


dr 


The  resulting  angles  of  deflection  are  rather  complicated  functions 

or  r and  b , so  that  it  would  be  difficult  to  consider  transitions  among 
m b 

the  elastic  trajectories.  Let  us  consider  several  approximate  trajectories 


which  would  facilitate  calculation.  We  observe  that  for  certain  rather  simple 
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trajectories  the  angle  of  deflection  is  uniquely  determined  by  the  ratio  r^/b  • 

For  r S b the  collision  is  repulsive,  for  r { b , it  is  attractive.  To 
m r r in  N 

investigate  the  possibility  of  approximating  the  true  trajectories  by  one  of  these 

simpler  cases  we  will  determine  for  trajectories  found  when  V OC 

i 

(Kepler  orbits).  V cC  y y (Maxwell  molecules)  and  for  a special  trajectory  of 
connected  line  segments.  Then  we  will  compare  the  functions  with 

values  for  the  Lennard  -Jones  potential.  For  the  Lennat  d - Jones  potential 
is  not  a single  valued  function  of  r^/b  . hut  comparative  numbers  can  be 
obtained  over  parts  of  the  range  by  considering  large  b for  attractive  collisions 
and  small  b for  repulsive  ones, 
a)  Kepler  trajectories 

A 

If  V - -p-  . the  equation  of  the  trajectory  is  easily  shown  to  be  the 


conic  section 


I 

r 


--A.  / 

v l *- 1 1 - e.  c-^/  6j 


(8) 


/7<rfc  x 

We  are  interested  in  the  repulsive  and  attractive  trajectories  which  are 
hyperbolae  described  by  equation  (8)  if  e ) 1 . If  k 0 , the  repulsive  case, 
the  center  of  force  turns  out  to  be  an  external  focus  of  the  hyperbola. 
y be  found  as  follows; 
from  equation  (8) 


u = 





(9) 
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from  equation  (7)  for  this  potential 

3 A ^ 

r. 


i - 


or 


*3 


- = O 


Jz_ 

"2 


_ J=_v  / JOzl 
/ 2. 


(10) 


solving  for  e , 


l^j  - / 


<2.  - 


^>r  ^ fc> 


(tr-  / 

are  the  directions  of  the  assymptotes, 

~ - — — — 

f e {^)  * / 


(ID 


(12) 


Since  the  angle  of  deflection  is  given  by 


we  have,  finally, 


<^02,  X 


(13) 


(14) 


For  attractive  collisions,  k ( 0 and  r b , the  center  of  force 

m 

•X 

becomes  an  internal  focus  of  the  hyperbola.  Letting  V ^ "p"",  and  proceeding 


as  ; bovq, 


from  equation  (8) 


I'T* 


M f h" 

Xk. ' ( I +-  e ) 
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and  from  equation  (7) 


/ b 


- o 


we  have 


/V  £ v TVm  ' ^ ^ 

c^r  - / 


e = 


2»o  as  before 


ft)'-/ 


CxT^X  ' / " 2. 


fei  - /_ 

(£Tw  J 


b)  Repulsive  potential  proportional  to  "py  (Maxwell  molecules) 

4c 

Let  V x "Ty  , then  the  equation  of  the  orbit  is 


Oo 

r 


0 - i [ XI  - r~- 


-y~ 

<±r 


I . - V' 


(15) 


where  E - M ^ . b is  the  collision  parameter  or  miss  distance,  and 
where  0 - 0 is  parallel  to  the  initial  asymptote. 


Let 


X - 


Xo  - bftY 


then 


X- 


- -L- 


H'-*v  'If. 


dX 


(16) 


d r 

The  apse  of  the  orbit,  where  -jrr  a 0 , is  given  by  the  real  positive  root  of 

O 0 


(17) 


Call  this  root  X / - ~p~  ■ Then 


X - 7r  - 
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The  integral  can  be  evaluated  in  terms  of  a complete  elliptic  integral  of  the 
first  kind. 

Let 


2 Xc 


_V<  *■  t- + '] 


(IS) 


»2*»2=  X o'  \ n- 


then 


X = 7 7 - 2xr 


(19) 


This  is  a standard  form  of  the  incomplete  elliptic  integral  of  the  first  kind*, 

that 

2 


X - 77  (h-ir-)'A  F(£,pJ 


where 


F, 


<p 


(X.  b 


Substituting  from  (17)  and  (18) 


X-  7T  - 


IX, 


v ^ - x ' 


(20) 


(21) 


Jahnke  and  Emde  Tables  of  Functions,  Dover,  1945. 
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whe  r e 


ir 

U , - f 

-j  — 


£ w )y 


the  complete  elliptic  integral  of  the  first  kind, 
c)  Attractive  potential  proportional  to  ' JTTj 


-Jk. 

r v 


Let  V ; ~^Z  < then  proceeding  as  for  Maswell  molecules. 


OO 

r t 


6 =-  - ±- 

J l b~  5 bv 


d r 


(22) 


Let 


A - 


Jl 
E b‘ 


then 


Letting 


gi  vet 


"O 


G = 


^iL 


oc'"  - X £ / ~ V~-  y*J 
X = a [_/-*-  V / - v/gj 


X - 7T  - 2. 


oo 


j 1 ITp^Xy^r) 


U3) 


where  y ^ is  the  largest  positive  real  root 


of 


y4  - y2  ♦ A = 0 


t.e.  y^  - J Again  we  have  a standard  form  of  the  complete  elliptic 
integral,  giving,  finally: 


X =-  11  ~ t Ht'tW) 


(24) 


W1S-ONR-  5 
9 October  1953 

14 


<1  • Str  e’ght  1 iu  segments 

Give;,  r -•  r-ci  b , wc  t art  sp pc  i f y ni  approximate  trajectory  such  that 
'he  angle  of  det’.cciion  d<  pertd'  oniv  on  r /b  by  proceeding  as  follows: 

For  repulsive  toll  s.ons  • ^ i>)  * he  molecule  continues  on  its 

o r ig.n  a ! M » a ght  1 .ne  path  uni  • 1 it  re  n hc->  a distant  e r ^ f rom  the  origin. 

\t  th  po.  • the  mole.  .U  - dt  i .<•  ted  e though  it  were  rebounding  from  a 
rigid  sphere  See  Fig.  1 ; . 


F.g  l 


Simple  geometry  gives  the  angle  of  drfle<tion 

co"  X - 2 (~j~-  / (25) 

For  attractive  collisions  <r  / b>  the  molecule  continues  on  the 

m ^ 

straight  line  unt  1 it  reac  hes  a distance  b from  the  origin.  At  this  point 
it  is  deflected  through  an  at  ute  angle  such  that  the  new  segment  passes  a 
d’  st  am  e r t rom  t be-  origin.  I he  remainder  of  the  trajectory  is  symmetric 
See  Fig.  L 


♦ 

I 

r.. 


m 


> 
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In  this  case 

cos  X - 2.  " / (26) 

In  the  following  graph  the  functions  cob  Xl  r^A. J are  plotted  against 

rm/b  for  the  various  trajectories  considered.  For  comparison,  some  Lannard- 
Jones  values  are  given  for  several  values  of  the  reduced  relative  kinetic  energy 

£ 

From  the  curves  of  Fig.  3,  it  is  clear  that  the  Kepler  orbits  are  the 
poorest  approximation  to  the  trajectories  produced  by  the  more  sharply  varying 
interanolecular  potentials. 


The  Molecular  Theor y of  Gases  and  Liquids  by  Hirschfe lde r,  Curtiss,  and 
Bird  (John  Wiley  It  Sons,  New  York,  1954)  (Appendix,  Ch.  6 with  Ellen  L. 
Spotz) 


5 
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Choice  between  the  line  segments  and  the  1 / r orbits  would  depend 

on  ease  of  manipulation  and  success  in  giving  the  desired  cross  sections  For 

repulsive  encounters  the  connected  line  segments  give  remarkably  good 

answers.  Since  these  line  segment  orbits  would  be  expected  to  give  the  largest 

deflection  of  any  purely  repulsive  trajectory  characterized  by  a certain  r , it 

m 

must  hr-  *!^t  the  attractive  part  of  the  Lenna  rd -Jones  potential  tends  to  compensate 
for  the  softness  of  the  repulsive  part,  resulting  in  a nearly  hard  sphere  deflection. 

In  attractive  encounters  the  line  segments  give  good  answers  for  r^ 
slightly  less  than  b , but  they  do  not  reflect  the  orbiting  tendency  as  r^  is 
decreased  The  line  segments,  for  exa  ..pie  give  full  180°  deflection  only  as 
approaches  zero. 

To  complete  the  specification  of  any  of  these  approximate  trajectories, 
we  need  only  assign  a speed  to  the  particle  at  each  point.  Clearly  this  cannot, 
in  general,  be  done  so  as  to  consec  ve  both  energy  and  angular  momentum. 

Since  we  plan  eventually  to  use  the  trajectories  in  a time  dependent  perturbation 
calculation  of  transition  probabilities,  we  arc  most  concei  led  with  the  time 
• pent  in  each  region  of  the  potential  Now,  the  approximate  trajectories  have 
been  adjusted  so  as  to  traverse  the  same  region  of  the  potential  as  the  true 
trajectory.  T-w,  f 4>°°  > 80  lt  would  appear  reasonable  to  adjust  the  speed  at 

each  point  so  as  to  keep  the  total  energy  of  the  system  constant.  Angular 
momentum  conservation  would  adjust  the  speed  according  to  the  direction  of 
motion  at  the  point  in  question,  a feature  of  little  interest  to  us  as  long  as  the 
total  deflection  of  the  path  is  nearly  correct. 
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It  is  ex  cauraging  to  note  that  at  the  d'  stanc  e of  closest  approach  our 
attractive  trajectories  coincide  wiMi  a point  c a true  trajectory  with  the 
correct  energy  and  angular  mementum  parameter  si  -It  is  from  this  region 
th&t  we  expect  the  greatest  contribution  to  transition  probabilities. 


C.  THE  INELASTIC  TRANSFER 


Consistent  with  our  perturbation  picture  wewill  w int  to  consider 
transitions  between  elastic  trajectories.  Before  investigating  such  transitions 
in  terms  of  any  of  the  approximate  trajectories  considered  in  the  previous 
section,  however,  let  us  look  at  the  limitations  imposed  by  the  conservation 
laws  of  classical  mechanics  on  any  tiajectory  at  the  moment  of  inelastic  transfer. 

The  dynamic  state  ol  an  encounter  between  linear  rotating  rigid 
molecules  in  conveniently  specified,  when  the  molecules  are  far  separated,  by 
the  eight  variables  b,  g,  ti  a , and  , where  we  have  chosen  a coordinate 

system  such  that  the  motion  lies  in  the  xy  plane.  / a and  ^b  are  the  angular 
momenta  due  to  the  rotation  of  molecules  a and  b . The  angular  momentum 
due  to  translational  motion  relative  to  the  center  of  mass  is  Jl- 7-  — tf  ^ b 
The  internal  energies  6 a and  £b  and  the  translational  energy  6 rare  easily 
determined  in  terms  of  these  variables. 

r^et  us  suppose  that  at  some  point  along  a trajectory  determined  by  these 
initial  conditions  and  the  potential  function  there  occurs  an  instantaneous 

change  in  the  internal  dynamic  state  such  that  a £'■(  and  Jb  -»^£b  with  the 
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new  potential  V , ^ Applying  the  conservation  laws  of  energy  and  angular 
mementum  we  have 

E = +6  b €-Krj  ♦ V’y.ri  — + £*  +■  €r(rJ  a-  V-r) 

(27) 

L ~ ^4.  + f A9r  = 7-  JLt 

To  interpret  the  consequences  of  these  equations,  consider  figure  4 


Here  we  picture  the  situation  in  terms  of  the  equivalent  one  particle  problem 
as  usual.  RS  , lying  in  the  xy  plane,  is  the  segment  of  the  trajectory  just 
before  the  transition.  The  transition  occurs  at  S , so  ST  is  the  trajectory 
immediately  after  transition.  In  general  ST  does  not  lie  in  the  xy  plane, 
but  in  a new  plane  determined  by  ST  and  point  0 ; we  designate  this  as 
plane  F The  xy  plane  and  F intersect  along  OS  . The  molecules  are 
separated  a distance  r at  the  instant  of  transition. 

Now,  the  total  angular  mementum  is  conserved,  so  in  particular  the 


! 
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component  along  OS  must  bo  icnsl.ii  t Thus: 

^ r i < U * c r j _ ^ ' 1 <■»  ^ Jcs 

, / 

But  and  tr  are  clearly  pcrpendu  ui.ir  to  OS  so  that 


. • i ~ ' £»/  -t. 

h 


' ( ti  ’ -ks-J  j 


/—  1 A/*- 


- y<JU<>x 


*1 


(28) 


or  where  ^ is  the  total  internal  angular  momentum,  and  f | is  the  angle  in  the 
xy  plane  as  shown  in  Fig  3 This  result  might  be  looked  upon  as  a restriction 
on  the  possible  redistribution  o(  x and  > lomponents  of  the  internal  angular 
mementum  for  a given  point  on  the  trajectory  Alternatively,  and  more  con- 
venient for  our  purposes,  it  specifies  the  point  or  points  (if  any)  on  the 
trajectory  at  which  a spec  lfied  transition  could  take  place.  Let  us  see  how  this 
restriction  might  operate  in  prac  tice  for  any  purely  attractive  or  repulsive 
trajectories.  If  0 - 0 is  the  direction  of  the  initial  asymptote  of  a trajectory, 
as  in  Fig.  4,  then  for  a repulsive  c ollision  the  direction  of  the  final  asymptote 
Gi  will  lie  in  the  range  £ 4 Yv  i-  ' . Since  the  tangent  function  is  single  valued  on 
the  interval  ( Gy  a specified  transition  c an  lead  to  at  most  one  allowable 
solution  of  equation  Z8  For  attractive  collisions,  on  the  other  hand, 
so  that  there  will  certainly  be  one  solution  J\  on  the  interval  ( 0t  7TJ  and  it 
is  quite  possible  that  1 7!  will  also  be  acceptable.  Qualitatively,  however,  it 
can  already  be  observed  that  there  *vil!  be  a rather  limited  range  of  Yj  such 
that  the  two  points  on  the  curve  represented  by  'n  and  he/Y  will  lie  at  about 


the  same  distance  from  the  origin.  Since  wc  t xpoct  that  the  probability  of 
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transition  will  fall  off  rapidly  with  increasing  separation  of  the  molecules,  it 
will  frequently  happen  that  only  one  of  the  allowed  transition  points  Is 
physically  important. 

From  the  energy  equation; 


6r0b)  - x ^ 


/ 

- £ - 6a  - 6b'  - W>)  . 


(29) 


Solving  for  g/  wo  have  the  speed  of  the  particle  after  the  transition.  If  this 
result  is  imaginary,  we  conclude  that  the  transition  is  energetically  forbidden. 

It  ie  interesting  to  note  at  this  point  that  if  one  were  considering  this  problem 
in  terms  of  curves  on  a multidimensional  energy  surfaces,  this  last  result 
worn!*  simply  be  the  Franck-C-ondon  principle  restricting  transitions  between 
sueh -curve*. 

To  complete  our  analysis  of  the  transition  we  use  the  angular  momentum 

equation  again: 

Jr  - L -M.  -41 

The  direction  of  £ determines  the  plane  of  the  new  trajectory.  The  selection 
rule  on  has  assured  us  that  this  plane  will  pass  through  5b  (Fig.  4)  so 
this  last  item  of  information  gives  us  the  angle  between  plane  F and  the  ay 
plans.  From  the  magnitude  of  Jj  we  determine  the  new  collision  parameter 

' =.  1 tr'I..  (M) 

9' 

Apparently,  if  b > E no  satisfactory  trajectory  can  be  constructed,  so 
this  result  constitutes  a third  selection  rule  on  the  transition.  If  b <.  , 
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surprisingly  enough,  there  are  two  possible  directions  for  the  new  trajectory, 
coriespondlng  to  the  two  tangents  to  a circle  of  radius  b which  pass  through 
a point  in  the  plane  outside  of  the  circle.  Each  of  these  directions  is  consistent 
with  the  total  angular  me  men  turn  of  the  system.  The  dual  choice  may  also  be 
seen  directly  from  the  expression  for  angular  mo  men  turn,  L * r xjfc  • Working 
in  the  plane  of  r and  p , JL  * rp  sin  (p  , where  & is  the  angle  from  £ to  . 
Apparently  for  a given  L , x • end  £ , both  f and  7 r - f>  are  acceptable 
solutions. 

That  the  choice  is  really  unique,  however,  may  bo  argued  as  follows: 

Let  us  suppose  that  the  inelastic  transition  specified  by  -»  and 

A Qjl 

approaches  aero  in  such  a way  that  - at  defined  by  equation  (28)  remains 

constant,  thus  keeping  the  . ranaltion  at  the  same  point  on  the  trajectory.  In  the 


limit  of  no  transition,  the  equations  of  motion  say  that  the  particle  continues 
on  the  smooth  course  rather  than  taking  the  abrupt  change  in  direction  which 
would  also  ho  consistent  with  the  energy  and  angular  memo  it  turn  of  the  syetem. 
Thus,  if  the  particle  is  on  the  approaching  branch  of  an  clastic  trajectory,  it 
does  mot  Jump  suddenly  to  the  corresponding  point  on  the  receeding  branch. 

Now,  a finite  transition  at  the  apecifitd  point  on  the  trajectory  may  be  considered 
at  a oorioaof  small  transitions  with  constant  » Then,  by  continuity  of  the 

whole  process,  the  particle  finds  itself  on  that  point  of  tho  new  trajectory  which 
correopoads  to  the  starting  point  on  tho  eld  trajectory,  correspondence  being 


Interpreted  in  terms  of  tho  picture  in  which  one  trajectory  is  continuously 
transformed  into  the  other  by  a progressively  increased  magnitude  of  the 
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inelastic  transfer.  In  practice  ‘his  simply  means  that  if  the  particle  is  on  the 
approaching  (receedlng)  part  of  one  fraiectorv,  it  must  appear  on  the  approarhi 
(receeding)  part  of  the  new  trajectory. 

Thus,  in  terms  of  a proposed  transition,  we  have  determined  where  the 
transfer  could  take  place  and  what  the  speed  and  direction  of  motion  of  the 
particle  must  be  after  the  transition  1 n the  course  of  getting  this  information 

from  the  energy  and  angular  momentum  conservation  laws,  we  have  found 
three  rather  simple  criteria  or  "selection  rules"  for  determining  whether  the 
transition  could  take  place  at  all: 

!.  ^ ..  given  by  equation  28  must  lie  on  the  initial  trajectory. 

2.  $ a.  given  by  equation  29  must  be  real. 

3.  b as  given  by  equation  30  must  be  less  than  the  separation  of 

the  molecules  at  the  time  of  transition. 

The  second  and  third  of  these  criteria  simply  tell  whether  or  not  the 
transition  has  taken  place  in  a region  compatable  with  the  total  energy  and 
angular  momentum  of  the  system.  The  first  is  a restriction  on  the  redistn  ■ 
button  of  angular  momentum  due  to  the  fact  that  the  trajectory  must  be 
continuous. 

In  the  next  section  we  will  consider  the  application  of  these  results  to 
the  problem  of  compounding  our  approximate  elastic  trajectories  to  represent 


inelastic  encounters. 
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I 

I D APPROXIMATE  TRAJECTORIES  FOR  INELASTIC  COLLISIONS 

i 

! 

To  be  consistent  with  our  approximations  at  this  point,  we  must  seek 
«n  < ne  la  >t  :c  t r ajec. « o t y sue  h that  energy  i fl  conserved  at  every  point  or.  it  and  such 
that  the  angular  momentum  is  conserved  when  the  molecules  are  far  separated 
before  dr.d  alter  the  encounter.  Such  a trajectory  may  be  constructed  by  the 
lolloping  procedure; 

l Given  ^ J b b g and  V ( r ) we  have  the  approximate  elastic 
trajectory  r T(<31  . Specifying  the  transition  by  the  new  rotational  angular 
momenta  / ^ and  *'*■'  solve  for  an  angle  >j  as  in  the  preceeding  section 

If  the  lines  £ ^ or  & ~ >n  — ^ do  not  interest  T , the  transition  is 

I 

impossible.  If  one  or  both  ol  these  intersections  exist,  we  have  possible  points 
for  the  transition.  Suppos.  r T.  , is  an  intersection. 

2.  At  rQ  wo  investigate  equation  29,  the  energy  equation,  to  determine 

the  new  relative  speed  g'  . Price  energy  is  to  be  conserved  at  every  point, 

' e 

we  appiv  the  second  sc  le  i t lon^  i go  r oua  ly  and  call  the  transition  forbidden  if 
g^s  imaginary.  Suppose  g/  is  .in  afrcptable  solution  of  equation  29. 

3.  Stric'  application  ot  the  third  selection  rule  on  angular  momentum 
it  the  transmon  point  would  now  be  inconsistent  Sm  h a restriction  would 
demand  that  this  point  on  the  new  trajectory  r T ( 0 ) have  correct  energy 
and  angular  momentum,  'oo  stringent  1 condition  for  an  appro*  i mate  trajectory 
Though  certainly  not  unde  s i r able  in  itself  sue  h a u strn  turn  a » t hi  s point 


would  complele'y  determine  the*  rest  of  the  path  md  would  not  leive  us  tree  to 
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adjust  the  asymtotic  value  of  the  angular  momentum. 

Instead,  we  simply  use  the  direction  of  Sr  - L to  give  us 

the  orientation  of  the  plane  F m which  the  new  trajectory  r = T* { &' ) must 

lie.  Since  the  elastic  trajectories  lie  in  a plar.e,  this  will  also  be  acceptable 

asymptotically.  From  E , / L/  and  V^r)  we  can  determine  the  shape  of 

T^d')  completely,  lacking  only  its  orientation  in  the  plane  F . Apparently 

if  r/  is  greater  than  r , no  orientation  of  T/  will  provide  a continuous 
m o 

transition  between  T and  , so  that  the  third  selection  rule  appears  here 

in  a rather  denatured  form.  If  t/  is  less  than  r there  will  be  two 

m o 

orientations  of  T/  in  F which  will  intersect  T at  r , but  from  the 

o 

argument  in  the  preceeding  section  the  choice  must  be  made  so  that  the 
particle  continues  approaching  or  receeding  from  the  apse  of  T*  according 
to  its  motion  on  T . 

Thus,  we  have  a simple  and  unique  method  for  construction  of 
approximate  inelastic  trajectories. 


E.  THE  INELASTIC  DIFFERENTIAL  SCATTERING  CROSS  SECTIONS 

So  far  we  have  described  a method  for  constructing  approximate 
inelastic  trajectories  in  terms  of  a combination  of  classical  and  quantum 
mechanical  concepts.  We  have  assumed  a classically  exact  knowledge  of  all 
coordinates,  but  at  the  same  time  we  considered  a discontinuous  change  in  the 
internal  coordinates  - a distinctly  quantum  mechanical  idea.  We  now  wish  to 
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i-  ho  t b s r.:o  • in-r 1 ;.r  1 r>  •.<* ; vc  t hi  rji;  in'iiin  int  h <i. : > ■*(  ‘Mori’  g problem  a nd 

, " * 

r\r  < r » rn  1 1 r I h 0 . T i I s " ! > d • ! • ' < I t ^ , • ,j  t <fj) 

: t i «o»m*  o!  * h > h l.'hionix-.  k ‘ > - d«>>  I ■ > p«  *1  ■>  theory  ot  the  »r  -tn^port  phenomena, 

» * ( ry  v,/  may  bo  donned  is  t o j i oa  s Lo*  a,  tr.d  >y  be  the  polar  and 

a r mo«  hal  angle  s of  the  i nil  • e 1 1 1 • --e  vc  lot  :t  v veco-  me  a sured  from  the 

m.t  •»!  re.attve  velor:*'-  vet  'o-  ' let  and  j « gmfv  al!  of  the  necessary  quantum 

numbf  r s to  specify  t Ke  :ri  » ; -» 1 internal  quint  um  states  o!  mo  le  c ules  a and  b' 

and  let  k and  1 represen-  iheHr  quantum  numbers  after  the  encounter.  Then 
A f 

! /u  cc  cued  tpi*  the  prob.ab'l'ty  »Kat  from  a unjfor  m be  am  (plane  wave! 

> d ' ' o 

of  molecules  :n  state  i incident  on  a molecule  in  state  j.  with  'elative  velocity 
g , a given  molecule  is  sca'tered  into  the  solid  angle  si  i*.  d&df  about 
while  the  interna!  quantum  states  go  to  k and  l during  the  encounter. 

If  J ^ and  4rff  expressed  in  tr-ms  of  the  polar  variables  q ^ , 

and  ^b  * 6 b ,*u’  Ru  *nt  um  m^5C  T 8 ’j  ^ are  sufficient  *o 

spec  >fy  the  c lassie  al  var  : able  h ^ ^ £ 4 rb  9 b,  y*  fj\  x b and  9 £ . 

Thus  tbe  definition  of  the  scattering  < ross  section  can  oe  reformulated  in  terms 
of  classical  variables  as  follows;  in  the  encounter  specified  by  g , A Q ^ 

, and  £b  . what  ta  the  probability  of  scattering  into  sine*  dd.  d V’  while  the 
internal  coerrdmates  go  to  0 ' , fj  ‘ , C.  . and  £ ’ 

In  order  to  specify  cx  and  -f*  , f’om  our  classical  derivation,  however, 
we  need  information  about  the  extra  variables  (p  , and  < p ^ or 

their  equivalent.  If  we  can  get  the  probability  of  the  desired  scattering  event 
by  assigning  probabilities  to  these  extra  va  • tables  according  to  the  quantum 
mechanical  correspondence  principle,  then  average  over  the  extta  variables. 


I 
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we  ought  to  have  an  approximation  of  the  differential  cross  sections,  valid  in 
the  correspondence  limit. 

For  a uniform  beam  of  incident  particles  (the  corresponden  e Unit  of 
a plane  wave)  the  relative  probability  of  a collision  parameter  in  the  range  b , 
db  is  just  77b  db  Similarly  in  the  correspondence  limit  the  azimuthal 
angle  of  the  angular  momentum  vector  of  a rigid  rotor  is  completely  uncertain, 
so  that  the  probability  density  function  for  the  unprimed  s is  just 
For  the  primed  ^ s , however,  the  situation  is  a little  different.  If  the  transi- 
tion occurs  at  a point  r on  the  trajectory  all  values  of  £p  a and  cp ^ are 

equally  likely  provided  they  are  consistent  with  the  first  selection  rule 
AjP* 

Apparently  this  problem  can  be  handled  if  we  change  variables 
from  (p ^ and  to  the  arc  length  s , according  to  the  first  selection  rule, 

and  a variable  z equal,  say,  to  cp ^ Now,  if  we  can  find  a function  ~f~(u) 
such  that  f (s)  ds  is  the  probability  of  transition  in  the  interval  s , ds,  then 
we  know  that  all  values  of  z will  have  probability  consistent  with  all  and 

Cp'h  equally  likely  If  z = p ^ , for  example,  all  z will  be  equally  probable 

f / 

since  the  relation  between  cp  and  p^  is  satisfied  automatically  by  the 
choice  of  s . 

The  function  ~f(s)  ds  is  easily  gotten  from  time  dependent  perturbation 

ft  , 

theory.  Let  us  define  the  function  T(ij  , kl  , b . g|  t)  so  that  -f  dt  is  the 
probability  that  the  ij  , kl  tranasition  takes  place  at  t , dt  when  we  are  on 
the  elastic  trajectory  determined  by  b , g , and  ij  . Let  and 

be  the  wave  functions  of  molecules  a and  b when  they  are  in  the  i'th  and 
j'th  quantum  states  respectively,  and  let  V(  Ya  , , t)  be  the  interaction 
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potential  at  t;mc  t . Then 

■t 


A/ 


% 


- <_  \ r ^ t - p K # 

* ! tl  / , 


(31) 


uhfrf  the  probability  of  tr.uw'  .on  up  to  fine  t is 


' * C ; :'j 


Here 


ClL:(z:)  - 0 k t i | 

- i ki.  i j 


and 


v“(t) 


r - * , .* 

! fa  §cl  =/_rb 


Dropping  the  extra  indices  tor  a moment,  let  us  expand  d (t)  in 


a Taylor's  series  about  some  time  tj  / 

d j 
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fhe  re  4 (»)  is  just  the  function  we  are  seeking  Explicitly 


f 


( ij  >-kC,  a,  $\t)  ~ 


2_ 

fiv 


t-l 


J 


fe* 


(32) 


Here  we  assume  that  the  encounter  started  at  t = -oo  Since  the  arc  length 
s is  a more  convenient  variable  for  us  than  the  time,  we  write 


{ f i b J It]  d t - 4 ( 1 1 **  *,  to  <i  I t(v) 

<±  S 

- 4 ( <■  1 , M.6,  b,  £ I Sj  d S 


dS 


(33) 


where  — is  simply  the  reciprocal  of  the  velocity  at  the  point  in  question. 

Before  writing  down  a general  formulation,  let  us  introduce  one  new 
function  J*  (lj  . , b , g , <p  a > <p ^ . ? . s)  which  symbolically  contains 

the  last  two  selection  rules.  Thus,  for  a point  s on  the  trajectory  determined 
by  b , g , lj  , and  ^ is  zero  if  the  lj  , kl  transition  is  forbidden  according 

to  the  last  two  selection  rules  and  is  unity  if  it  is  allowed. 

We  are  now  in  a position  to  write  down  a formal  procedure  for  determin- 
ing  ^ ■ In  actual  calculations  a judicious  choice  of  variables  and 

sequence  of  calculation,  as  well  as  some  approximations,  will  doubtless  present 
themselves.  The  following  formulation,  therefore,  is  to  be  taken  only  as  a 
formal  presentation  of  the  procedure. 

The  event  of  interest,  a collision  with  definite  scattering  angles,  is 
specified  classically  by  the  fourteen  variables  b , g , ij  , kl  , <p  , , 

(p a For  a given  g and  ij  , kl  we  wish  to  find  the  probability 
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dlstr lbution  function  of  the  variables  oC  and  Ip  . Let  us  use  the  first 
selection  rule  to  introduce  the  arc  length  ■ as  a new  variable  in  place  of 

‘Pi,  • 

For  these  variables  we  know  all  of  the  necessary  probability  functions. 


Xi 

Let  I ( §/  b,  fa.  P v . <fH,  S)  di  d<?e  d<pa  d£ 


>e  the  probability  of  an  inelastic  collision  in  which  the  variables  b , qP 

1 a 

( p ^ , (p'^  , and  s lie  in  the  specified  ranges.  From  the  preceeding 


discussion 

Kj  (%l  b 

— C-  ~1~(  ‘ ;«  5CfcL-*/.  b,£.  ¥*.¥•>,<&',  s) 

Now  for  a given  g . lj  , kl  let  us  introduce  new  variables 


(34) 


oC  = Afo.ja.  b f <P«'.S) 

£/  b,<P*  ,fo,  9*',s) 

J^(  = Xl(M.  *«,§/  b,  <j*\  SJ 

Xa.  = Xa.( s) 


X3  = X3(  <Pa,S) 
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Merc  X | X and  X • r,‘  proper! » i hotter:  ( ombin.it  ions  of  the  original 

vimbif  s vi  hi  i h tn  a r <*  t In'  tr  ui  s in » in  • * ion  « or.  vr  merit  and  non  singular  Very 
likely  the  \ , *:1!  brs1  be  > ho  sen  ,n  i s.inple  intm.fr  such  vs  X b 

Xz  * *"'i  X,  •?, 

Int  r odui  : n g hese  Me*-  . >r>ab!e->  in  tnsf  probability  function 

S;  it  d r1--  " to  c Pc*  aS 


rS  * 

V > I ^ 'r-'-  'x  c 


- - ? ’ - *•«>  <»> 


- ry«e 

~~  1 i • ( & I x.  X- X*-  Xi)  ck*.  d ip  d X ci  X v.  d Xj 


if  we  no*  integrate  over  the  appropriate  ranges  of  X.  \ • X 2 ’ dnc*  ^3  WC  ^'ave 


// 

p.r/ 

I C S 


; s ^ a W 


r r r 

ij! 

X Xv  t| 


I l ) ' 


X V ,‘i)  rtX. 


(37) 


No,  P"'§*  y)  d<  a / '3  jus*  the  probability  of  scattering  into  <*,  d°(  and  1^; 
so  that 


r X* 

n ■ >5  d d.  d 1/ 


— ‘ i ' ^ ^ xJ-a^ri  o(  cfd.  d 


or 


// 


_./</  L_ 


Thus  we  finally  have  the  de  s red  result: 


*4 


— t t ^ ^ . Vj  - 


77  oC 


(38) 

J|jb  dXl  dX' dXj 


